Comparing surface-soil moisture from the SMOS mission and the ORCHIDEE land-surface model over the Iberian Peninsula

The aim of this study is to compare the surface soil moisture (SSM) retrieved from ESA's Soil Moisture and Ocean Salinity mission (SMOS) with the output of the ORCHIDEE (ORganising Carbon and Hydrology In Dynamic EcosystEm) land surface model forced with two distinct atmospheric data sets for the period 2010 to 2012. The comparison methodology is first established over the REMEDHUS (Red de Estaciones de MEDición de la Humedad def Suelo) soil moisture measurement network, a 30 by 40. km catchment located in the central part of the Duero basin, then extended to the whole Iberian Peninsula (IP). The temporal correlation between the in-situ, remotely sensed and modelled SSM are satisfactory (r. >. 0.8). The correlation between remotely sensed and modelled SSM also holds when computed over the IP. Still, by using spectral analysis techniques, important disagreements in the effective inertia of the corresponding moisture reservoir are found. This is reflected in the spatial correlation over the IP between SMOS and ORCHIDEE SSM estimates, which is poor (¿. ~. 0.3). A single value decomposition (SVD) analysis of rainfall and SSM shows that the co-varying patterns of these variables are in reasonable agreement between both products. Moreover the first three SVD soil moisture patterns explain over 80% of the SSM variance simulated by the model while the explained fraction is only 52% of the remotely sensed values. These results suggest that the rainfall-driven soil moisture variability may not account for the poor spatial correlation between SMOS and ORCHIDEE products.Peer ReviewedPostprint (published version

Crop Yield Estimation and Interpretability With Gaussian Processes

This work introduces the use of Gaussian processes (GPs) for the estimation and understanding of crop development and yield using multisensor satellite observations and meteo- rological data. The proposed methodology combines synergistic information on canopy greenness, biomass, soil, and plant water content from optical and microwave sensors with the atmospheric variables typically measured at meteorological stations. A com- posite covariance is used in the GP model to account for varying scales, nonstationary, and nonlinear processes. The GP model reports noticeable gains in terms of accuracy with respect to other machine learning approaches for the estimation of corn, wheat, and soybean yields consistently for four years of data across continental U.S. (CONUS). Sparse GPs allow obtaining fast and compact solutions up to a limit, where heavy sparsity compromises the credibility of confidence intervals. We further study the GP interpretability by sensitivity analysis, which reveals that remote sensing parameters accounting for soil moisture and greenness mainly drive the model predictions. GPs finally allow us to identify climate extremes and anomalies impacting crop productivity and their associated drivers

L-Band Vegetation optical depth and effective scattering albedo estimation from SMAP

Over land the vegetation canopy affects the microwave brightness temperature by emission, scattering and attenuation of surface soil emission. Attenuation, as represented by vegetation optical depth (VOD), is a potentially useful ecological indicator. The NASA Soil Moisture Active Passive (SMAP) mission carries significant potential for VOD estimates because of its radio frequency interference mitigation efforts and because the L-band signal penetrates deeper into the vegetation canopy than the higher frequency bands used for many previous VOD retrievals. In this study, we apply the multi-temporal dual-channel retrieval algorithm (MT-DCA) to derive global VOD, soil moisture, and effective scattering albedo estimates from SMAP Backus-Gilbert enhanced brightness temperatures posted on a 9 km grid and with three day revisit time. SMAP VOD values from the MT-DCA follow expected global distributions and are shown to be highly correlated with canopy height. They are also broadly similar in magnitude (though not always in seasonal amplitude) to European Space Agency Soil Moisture and Ocean Salinity (SMOS) VOD. The SMOS VOD values are based on angular brightness temperature information while the SMAP measurements are at a constant incidence angle, requiring an alternate approach to VOD retrieval presented in this study. Globally, albedo values tend to be high over regions with heterogeneous land cover types. The estimated effective scattering albedo values are generally higher than those used in previous soil moisture estimation algorithms and linked to biome classifications. MT-DCA retrievals of soil moisture show only small random differences with soil moisture retrievals from the Baseline SMAP algorithm, which uses a prior estimate of VOD based on land cover and optical data. However, significant biases exist between the two datasets. The soil moisture biases follow the pattern of differences between the MT-DCA retrieved and Baseline-assigned VOD values

Comparison of measured brightness temperatures from SMOS with modelled ones from ORCHIDEE and H-TESSEL over the Iberian Peninsula

L-band radiometry is considered to be one of the most suitable techniques to estimate surface soil moisture (SSM) by means of remote sensing. Brightness temperatures are key in this process, as they are the main input in the retrieval algorithm which yields SSM estimates. The work exposed compares brightness temperatures measured by the SMOS mission to two different sets of modelled ones, over the Iberian Peninsula from 2010 to 2012. The two modelled sets were estimated using a radiative transfer model and state variables from two land-surface models: (i) ORCHIDEE and (ii) H-TESSEL. The radiative transfer model used is the CMEM. Measured and modelled brightness temperatures show a good agreement in their temporal evolution, but their spatial structures are not consistent. An empirical orthogonal function analysis of the brightness temperature's error identifies a dominant structure over the south-west of the Iberian Peninsula which evolves during the year and is maximum in autumn and winter. Hypotheses concerning forcing-induced biases and assumptions made in the radiative transfer model are analysed to explain this inconsistency, but no candidate is found to be responsible for the weak spatial correlations at the moment. Further hypotheses are proposed and will be explored in a forthcoming paper. The analysis of spatial inconsistencies between modelled and measured TBs is important, as these can affect the estimation of geophysical variables and TB assimilation in operational models, as well as result in misleading validation studies

Impact of day/night time land surface temperature in soil moisture disaggregation algorithms

Since its launch in 2009, the ESA’s SMOS mission is providing global soil moisture (SM) maps at ~40 km, using the first L-band microwave radiometer on space. Its spatial resolution meets the needs of global applications, but prevents the use of the data in regional or local applications, which require higher spatial resolutions (~1-10 km). SM disaggregation algorithms based generally on the land surface temperature (LST) and vegetation indices have been developed to bridge this gap. This study analyzes the SM-LST relationship at a variety of LST acquisition times and its influence on SM disaggregation algorithms. Two years of in situ and satellite data over the central part of the river Duero basin and the Iberian Peninsula are used. In situ results show a strong anticorrelation of SM to daily maximum LST (R˜-0.5 to -0.8). This is confirmed with SMOS SM and MODIS LST Terra/Aqua at day time-overpasses (R˜-0.4 to -0.7). Better statistics are obtained when using MODIS LST day (R˜0.55 to 0.85; ubRMSD˜0.04 to 0.06 m3 /m3 ) than LST night (R˜0.45 to 0.80; ubRMSD˜0.04 to 0.07 m3 /m3 ) in the SM disaggregation. An averaged ensemble of day and night MODIS LST Terra/Aqua disaggregated SM estimates also leads to robust statistics (R˜0.55 to 0.85; ubRMSD˜0.04 to 0.07 m3 /m3 ) with a coverage improvement of ~10-20 %.Peer ReviewedPostprint (published version

Multi-temporal evaluation of soil moisture and land surface temperature dynamics using in situ and satellite observations

Soil moisture (SM) is an important component of the Earth’s surface water balance and by extension the energy balance, regulating the land surface temperature (LST) and evapotranspiration (ET). Nowadays, there are two missions dedicated to monitoring the Earth’s surface SM using L-band radiometers: ESA’s Soil Moisture and Ocean Salinity (SMOS) and NASA’s Soil Moisture Active Passive (SMAP). LST is remotely sensed using thermal infrared (TIR) sensors on-board satellites, such as NASA’s Terra/Aqua MODIS or ESA & EUMETSAT’s MSG SEVIRI. This study provides an assessment of SM and LST dynamics at daily and seasonal scales, using 4 years (2011–2014) of in situ and satellite observations over the central part of the river Duero basin in Spain. Specifically, the agreement of instantaneous SM with a variety of LST-derived parameters is analyzed to better understand the fundamental link of the SM–LST relationship through ET and thermal inertia. Ground-based SM and LST measurements from the REMEDHUS network are compared to SMOS SM and MODIS LST spaceborne observations. ET is obtained from the HidroMORE regional hydrological model. At the daily scale, a strong anticorrelation is observed between in situ SM and maximum LST (R ˜ -0.6 to -0.8), and between SMOS SM and MODIS LST Terra/Aqua day (R ˜ - 0.7). At the seasonal scale, results show a stronger anticorrelation in autumn, spring and summer (in situ R ˜ -0.5 to -0.7; satellite R ˜ -0.4 to -0.7) indicating SM–LST coupling, than in winter (in situ R ˜ +0.3; satellite R ˜ -0.3) indicating SM–LST decoupling. These different behaviors evidence changes from water-limited to energy-limited moisture flux across seasons, which are confirmed by the observed ET evolution. In water-limited periods, SM is extracted from the soil through ET until critical SM is reached. A method to estimate the soil critical SM is proposed. For REMEDHUS, the critical SM is estimated to be ~0.12 m3/m3 , stable over the study period and consistent between in situ and satellite observations. A better understanding of the SM–LST link could not only help improving the representation of LST in current hydrological and climate prediction models, but also refining SM retrieval or microwave-optical disaggregation algorithms, related to ET and vegetation status.Peer ReviewedPostprint (published version

Extended triple collocation: estimating errors and correlation coefficients with respect to an unknown target

Calibration and validation of geophysical measurement systems typically require knowledge of the true value of the target variable. However, the data considered to represent the true values often include their own measurement errors, biasing calibration, and validation results. Triple collocation (TC) can be used to estimate the root-mean-square-error (RMSE), using observations from three mutually independent, error-prone measurement systems. Here, we introduce Extended Triple Collocation (ETC): using exactly the same assumptions as TC, we derive an additional performance metric, the correlation coefficient of the measurement system with respect to the unknown target, rho(t,Xi). We demonstrate that rho(2)(t,Xi) is the scaled, unbiased signal-to-noise ratio and provides a complementary perspective compared to the RMSE. We apply it to three collocated wind data sets. Since ETC is as easy to implement as TC, requires no additional assumptions, and provides an extra performance metric, it may be of interest in a wide range of geophysical disciplines.Peer ReviewedPostprint (published version

Surface moisture and temperature trends anticipate drought conditions linked to wildfire activity in the Iberian Peninsula

In this study, drought conditions involving risk of fires are detected applying SMOS-derived soil moisture data and land surface temperature models. Moisture-temperature (SM-LST) patterns studied between 2010 and 2014 were linked to main fire regimes in the Iberian Peninsula. Most wildfires burned in warm and dry soils, but the analysis of pre-fire conditions differed among seasons. Absolute values of SM-LST were useful to detect prone- to-fire conditions during summer and early autumn. Complementarily, SM-LST anomalies were related to droughts and high fire activity in October 2011 and February-March 2012. These episodes were coincident with abnormally anticyclonic atmospheric conditions. Results show that combined trends of new soil moisture space-borne data and temperature models could enhance fire risk assessment capabilities. This contribution should be helpful to face the expected increase of wildfire activity derived from climate change.Peer ReviewedPostprint (published version

How many parameters can be maximally estimated from a set of measurements?

©2015 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.Remote sensing algorithms often invert multiple measurements simultaneously to retrieve a group of geophysical parameters. In order to create a robust retrieval algorithm, it is necessary to ensure that there are more unique measurements than parameters to be retrieved. If this is not the case, the inversion might have multiple solutions and be sensitive to noise. In this letter, we introduce a methodology to calculate the number of (possibly fractional)Peer ReviewedPostprint (published version

Relationship between vegetation microwave optical depth and cross-polarized backscatter from multiyear Aquarius observations

Soil moisture retrieval algorithms based on passive microwave remote sensing observations need to account for vegetation attenuation and emission, which is generally parameterized as vegetation optical depth (VOD). This multisensor study tests a new method to retrieve VOD from cross-polarized radar backscattering coefficients. Three years of Aquarius/SAC-D data were used to establish a relationship between the cross-polarized backscattering coefficient σHV and VOD derived from a multitemporal passive dual-channel algorithm (VODMT). The dependence of the correspondence is analyzed for different land use classes. There are no systematic differences in the slope for woody versus nonwoody vegetation, resulting in a strong correlation (80% explained-variance) and a global linear relationship when all classes are combined. The relationship is stable over the years of observations. The comparison of the Aquarius-derived VODMT to Soil Moisture and Ocean Salinity's multi-angular VOD estimates shows similar spatial patterns and temporal behavior, evident in high correlations. However, VODMT has considerably higher mean values, but lower dynamic range globally. Most of the differences can be attributed to differences in instrument sampling. The main result of this study, a relationship between backscatter and VOD, will permit high-resolution mapping of VOD with synthetic aperture radar measurements. These maps allow future studies of scaling and heterogeneity effects of vegetation on soil moisture retrieval at the coarser scales of land microwave radiometry. The study shows that VOD based on passive measurements and predicted by active measurements are comparable globally and that the breakdown by land cover classification does not affect the relationship appreciably